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SUMMARY 


The experimental investigation, conducted in 1970 under the 
auspices of the Army Air Mobility R&D Laboratory and the 
National Aeronautics and Space Administration, explored the 
mechanisms of retreating blade stall at one spanwise location 
on an eight-foot diameter model rotor blade. This test pro- 
gram provided blade airloads and motions data for a wide range 
of rotor operating conditions that could be used to evaluate 
current theoretical representations of the complex rotor en- 
vironment. As part of an additional contract awarded in 1971, 
a theoretical study based on this test data has been completed 
and is described in this report. The hybrid computer uniform 
downwash rotor simulation used in this study combines the 
storage and logic capabilities of digital equipment with the 
versatility of analog simulations to accommodate the complex 
elastic blade analysis and unsteady aerodynamics of the heli- 
copter rotor blades. 

The primary difference between the analog simulation and typi- 
cal digital formulations is in the number of radial integra- 
tions or updates completed in a rotor revolution. For the 
typical digital program, spanwise integrations are completed 
every 15° of azimuth, whereas, for the hybrid simulation, 
radial updates occur approximately every degree of azimuth. 

A comparison of experimental rotor performance data with 
hybrid simulation results at an advance ratio of .35 indicates 
that the theoretical formulation adequately predicts the rotor 
performance below stall but underpredicts the total lifting 
capability of the rotor beyond lift stall. For the advance 
ratio of .6, the theoretical predicted performance levels 
match the test data over the range of available thrusts. At 
the lower advance ratios, however, the limitations of uniform 
downwash imposed by the present hybrid formulation signifi- 
cantly affect the performance prediction capability. 

To further evaluate current theoretical representations of the 
complex rotor environment, a study of pertinent parameters as 
a function of azimuth and span was begun. Out of this study, 
three significant points arose. 

1. In the stalled regime, the hybrid prediction of blade 
flapping represents a significant improvement over the 
digital formulations examined. 
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. The present representations of the unsteady aerodynamic 
lift and moment in the stall regime are inadequate. 

3. In spite of deficiencies in the unsteady aerodynamic rep- 
resentation, the trend of the predicted growth in blade 
root torsion agrees well with test results. 

The improvement in the predicted flapping correlation over 
previous digital formulations is significant in that the 
proper location of the rotor in space is a prerequisite for 
adequate prediction of rotor airloads data. The improved pre- 
diction capability of the hybrid simulation is attributed to 
the increased number of radial updates inherent in the analog 
system that could significantly affect the converged solution 
to the equations of motion. A study has been suggested to 
determine the effect of update time on digital analyses. 

The deficiencies in the prediction of dynamic lift and moment 
levels possibly results from the manner in which two-dimen- 
sional oscillating airfoil data is used in the formulation of 
the dynamic stall delay. Most oscillating airfoil tests 
achieve representative high rates of change of oc by high fre- 
quency, low amplitude sinusoidal pitching motions about moder- 
ate values of oc . Data presented in this report suggests 
that more representative large amplitude pitch oscillations 
would greatly affect the value of the dynamic moment coeffi- 
cient and improve the prediction of blade element loads in the 
stall region. 

To evaluate the influence of scale effects in the model rotor 
test results, two separate decks of two-dimensional airfoil 
data representing model scale and full scale Reynolds number 
environments were input into the hybrid simulation. The hy- 
brid simulation, using full scale airfoil data, accurately 
predicted model blade loads. Using the model scale Reynolds 
number airfoil data, the hybrid formulation does not properly 
predict model blade loads. This suggests that the effective 
Reynolds number of the model rotor, as it affects the blade 
lift, is more nearly that of a full scale rotor. 

A preliminary study aimed at developing dynamic design cri- 
teria for airfoils indicates that modifying airfoil stall 
delay characteristics for negative values of oc can produce 
overall improvements in rotor performance and airloads. Since 
differences in the stall characteristics for negative c< and 
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the related damping characteristics can be changed by proper 
choice of airfoil, the relative benefits of the high static 
lift capability of leading edge stall airfoils can be tr jaded 
off against the benefits of the low negative damping. of trail- 
ing edge stall airfoils. Airfoil studies should continue with 
the aim of tailoring airfoils to produce improved stall re- 
covery behavior without seriously limiting dynamic lift capa- 
bility. 
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INTRODUCTION 


Helicopter forward speed and rotor design requirements have 
been dictated by limitations imposed by the effects of rotor 
blade stall. Thus, a better understanding of the stall phe- 
nomenon and a prediction capability representative of the com- 
plex rotor environment is imperative to the design of next 
generation helicopters. For many years, computer simulations 
of the helicopter rotor have not accurately predicted rotor 
performance or airloads above stall since they lacked proper 
dynamic stall characteristics and a proper representation of 
the three-dimensional environment of the rotor. 

In 1969, Harris, et al, empirically combined the airload char- 
acteristics associated with yawed flow and of a two-dimensional 
airfoil oscillating in pitch, and generated an advanced rotor . 
analysis method. These effects, when integrated over the rotor 
disc improved the rotor performance prediction capability, but 
no assessment of the capability of predicting local blade ele- 
ment airloads could be made without suitable test data for 
comparison. To provide such data. The Boeing Company, under 
contract to NASA, investigated the effects of stall on a rotor 
blade element in a three-dimensional rotating environment. 

The resulting model rotor test program provided blade element 
airloads and local boundary layer flow characteristics at the 
three-quarter blade radius position for a wide range of rotor 
operating conditions. A description of this test program and 
a summary of test results is given in Reference 1. The empha- 
sis of this referenced report was primarily on test results 
and provided only a limited study of theory/test correlation. 
Therefore, an additional contract was awarded in 1971 to pro- 
vide for final reduction and presentation of test data, and 
for a more detailed theoretical study. An informal data 
report containing final test data has been completed and may 
be obtained from the program technical monitor, W.J. McCroskey. 
The theoretical study, completed on the Vertol Hybrid Computer 
System, and the theory test comparison are the subjects of 
this report. 
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HYBRID STUDY 


In the first quarter of 1972, a theoretical study was under- 
taken to investigate the blade loads and motions associated 
with a three-bladed articulated rotor operating in the blade 
stall regime. Data points were chosen to match test condi- 
tions recorded in the 1970 eight foot pressure blade wind 
tunnel test program BVW 054. 

The hybrid simulation used in this theoretical study is a com- 
bination of analog and digital capabilities in that the airfoil 
characteristics and blade structural mode shapes are supplied 
to the analog equipment by a special purpose digital computer 
acting as a function generator. In addition, the logic and 
switching circuits required for the multiple integrations of 
this problem are provided by digital components. A description 
of the hybrid simulation is given in Figure 1. 

The formulation of the analog analysis is one of the uniform 
downwash with three blade flapping modes and two blade torsion 
modes. The representations ofunsteady aerodynamics and sweep 
effects are essentially the same as used in Vertol digital 
programs and as developed in Reference 2. The aerodynamic 
loads, and flapping and torsional motions resulting from these 
loads, are computed as they act radially along the blade. 

These blade loads are transferred to the appropriate coordinate 
system and doubly integrated to produce total rotor forces. A 
detailed description of the analog formulation and its capa- 
bilities is given in Reference 3. 

In the analog representation, the theoretical equations are 
transformed so that time is used to represent the azimuthal 
rotation of the blade as well as the radial sweep. The azi- 
muthal integration, however, proceeds much more slowly them the 
radial integration. One radial integration is completed within 
one degree of azimuth rotation. The airloads and motions com- 
puted during this radial sweep are held and used in the azi- 
muthal integrations until the following radial sweep. The 
major difference, therefore, between this hybrid simulation and 
standard digital programs is in the number of radial updates 
per rotor revolution. For typical digital programs, the span- 
wise integrations are completed every 15° of azimuth, and/ for 
the hybrid simulation, radial updates occur approximately 
every degree of azimuth. The significance of the increased 
number of updates is evident when considering the formulation 
with higher harmonics. Since many cycles are completed for 
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FIGURE 1 


HYBRID ROTOR AIRLOADS AND PERFORMANCE ANALYSIS 


MODEL 




1. Visual output of results 

2. Short run time 

3. On-line modified formulation 

4. Simplified parametric studies 

5. Adaptable to a variety of rotor types - Articulated, 

teetering, and rigid rotors 

6. Three flap modes, two torsion modes 

7. Provision for coupled modes 

8. Contains yawed flow and unsteady aerodynamics 
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each rotor revolution, increments smaller than the standard 
digital 15° increment are required to provide adequate repre- 
sentation of these higher harmonics. Such an increase in the 
number of radial updates for the digital computer would greatly 
affect the computation time and thereby the cost. 
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THEORY/ TEST COMPARISON 


Before a detailed study of the blade element airloads was be- 
gun, an assessment of the hybrid simulation capability of pre- 
dicting overall rotor performance was made. The data points 
for this study were chosen to match test conditions recorded 
during the 1970 eight foot pressure blade wind tunnel test, 

BVWT 054. These data points included both stalled and un- 
stalled conditions at advance ratios of .15, .35 and .6. In- 
put parameters for the hybrid study were based on model rotor 
structural properties and test program control positions. At 
a given advance ratio, collective pitch was set to match a given 
test thrust and then shaft angle was varied to produce a thrust 
sweep. A comparison of theory and test thrust sweeps is given 
in Figure 2 for an advance ratio of .35. The theory indicates 
good agreement with test up to lift stall, but underpredicts 
total rotor lift beyond lift stall. This difference in rotor 
lift level above lift stall is evident in the curves of normal 
force coefficient versus azimuth which are discussed below. 

The additional data shown in Figure 2 represents results from 
a similar rotor tested under the same conditions in the 
University of Maryland Wind Tunnel except for the difference in 
collective pitch. The similarity in the trends of this test 
data and the theory leads to more confidence in the capability 
of predicting rotor performance. The comparison of theoreti- 
cal and test propulsive force data in Figure 3 indicates fairly 
good agreement considering the magnitude of the hub tare cor- 
rection applied to the test ( data and the trend similar to the 
UMWT 516 test data. Curves of BVWT 054 test data are shown 
with and without hub tares to emphasize the size of the hub; 
tare correction. 

Similar theory/ test comparisons of rotor lift versus shaft 
angle at advance ratios of .6 and .15 are shown in Figure 4. 

At M - .6, the theory/test correlation is very good for the two 
different collective curves. For the lower advance ratio, a 
significant difference in the rotor lift curve slopes of the 
theory and test data is noted. This difference may be related 
to the significant nonuniformity in induced velocity along the 
longitudinal axis of the rotor as described by Harris in 
Reference 4. 

In this reference, the existence of the wake nonuniformity is 
indicated by large values of lateral flapping produced at low 
advance ratios. The effect of this nonuniform downwash on 
rotor performance is minimized at large negative control axis 
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angles (small wake skew angles). However, as the control axis 
angles are increased (increasing wake skew angles) , the non- 
uniform downwash tends to reduce the lift on the aft portion Of 
the rotor disc and produce an overall reduction in the total 
rotor lift curve slope. The comparison of the uniform down- 
wash theoretical solution with the nonuniform downwash test 
data shown in Figure 4 clearly illustrates such a change: in 
overall rotor lift curve slope. 

Encouraged by the reasonable agreement between the hybrid 
simulation and test oyerall rotor performance, at advance 
ratios of .35 and .6, a study of predicted blade root flapping 
was made to further evaluate the applicability, of the hybrid 
simulation. In typical digital formulations of the rotor en- 
vironment, the inability to predict blade flapping has been a 
serious shortcoming. Without proper representation of flap- 
ping, a theoretical formulation does not properly locate the 
rotor in space and, therefore, cannot provide the required 
velocity or angle of attack environment. Any attempted study 
of the resulting blade loads could prove to be misleading. 

Figure 5 provides a comparison of the hybrid root flapping pre- 
diction with the results of two current digital formulations 
(Vertol Programs B-67 and C-60) at an advance ratio of .35. 

The B-67 digital program is a rotor airload and performance 
analysis representing the blade by a series of blade segments 
producing a multi- filament near field, mid field, and far field 
vortex pattern and nonuniform induced inflow. The modal analy- 
sis utilizes two flapping modes and one elastic torsion mode, 
without flap-pitch coupling. Further description of this pro- 
gram is provided in Reference 5. The Vertol C-60 digital pro- 
gram employs a lumped mass model of the rotor providing results 
up to ten harmonics. The flap and pitch motions of the blade 
are coupled and the hub and lower controls are included in the 
analysis. The present wake representation includes only tip 
and root vortices. The c-60 data presented in this report was 
based on a uniform downwash program option described in Refer- 
ence 6. The unsteady aerodynamics formulations of the Hybrid, 
B-67* and C-60 simulation are essentially the same. The com- 
parison presented in Figure 5 clearly points out the short- 
comings of present digital formulation in positioning the rotor 
in space. The hybrid results, although still short of complete 
agreement, provide the best correlation to date with test 
flapping data . 
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Similar improvements in flapping correlation were noted for all 
thrust levels at the advance ratios of .35 and .6. For the 
lower advance ratio , M = .15, however, the hybrid solution 
produces no significant improvement. The differences in the 
predicted and test recorded thrust levels discussed above, re - 
lating to the limitations of uniform downwash, significantly 
affect the prediction capability at this advance ratio. It 
should be pointed out here that no systematic investigation was 
undertaken at the time of this study to improve the hybrid . 
flapping correlation further; however, such a study could prove 
to be beneficial. 

Since the aerodynamic formulations of the B-67, C-60 and hybrid 
simulations are very similar, the differences in flapping re- : 
suits could be attributed to the manner in which the dynamic 
characteristics are considered. The present digital formula- 
tions are limited in the number of radial calculations per 
rotor revolution to fifteen because of the run time and cost 
imposed by smaller increments. The increased number of radial 
updates inherent in the hybrid simulation (360 per revolution) 
could significantly affect the converged solution to the equa- 
tions of motion and the proper representation of higher elastic 
modes. It is suggested that a study be conducted on the ef- 
fects of updating time on blade loads and motions and particu- 
larly on root flapping. This could be accomplished by changing 
the number of radial computations performed by the B-67 digital 
analysis in a given rotor revolution. 

A comparison of predicted blade element airloads and motions 
with test results obtained from the pressure blade program 
provides an informative overview of the theoretical capability 
of this hybrid program. The curves of Figures 6 through 9 pro- 
vide theory/test comparison of blade element normal force co- 
efficient, moment coefficient and blade pitch at the three- 
quarter radius station. The predicted levels of normal force 
and moment coefficient are significantly underpredicted on the 
retreating portion of the rotor disc. There are indications 
that similar shortcomings exist in other theoretical, analyses. 
Reference 7. The theoretical formulation of the unsteady 
aerodynamic nature of the flow relies on a representation of 
the dynamic lift and moment based on oscillating airfoil test 
data, McCroskey points out in Reference 8 that the present 
computations of lift and moment overshoot are based on oscilla- 
ting airfoil data which are not representative of the azimuthal 
variation of oc exhibited by the rotor. Most oscillating air- 
foil studies have been conducted using high frequency, low 
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amplitude sinusoidal motions about small mean values of ot . 
The angle of attack variation with azimuth for a rotor more 
nearly resembles a large amplitude sine wave at the fundamen- 
tal rotor frequency. An investigation of available oscillating 
airfoil data including oscillations up to + 10° indicates that 
the value of dynamic moment is influenced by the amplitude of 
oscillation. Figure 10 presents oscillating airfoil data from 
Reference 9 at mean angles of 10° and 15° , ' illustrating a 
definite increase in minimum moment coefficient with increases 
in the amplitude of oscillation. Further study into the ef- 
fects of the oscillating airfoil pitch amplitude on the theo- 
retical unsteady aerodynamic y function representation and the 
resulting effects on dynamic lift and moment should be made to 
improve existing theoretical models. 

It is generally recognized that the moment stall indicated in 
Figure 7 has a significant effect on rotor control loads and 
in particular on the blade torsion loads. The influence of 
the large negative pitching moments in stall can be seen in 
Figure 9 since the elastic twist is linearly proportional to 
the blade torsion at the root of the .blade. The large pitch- 
ing moment generated at 240® is reflected as a nose down 
pitching motion of approximately 10°. With the following re- 
duction in -Cm, the blade begins a rapid pitch up motion 
followed by several oscillatory cycles commonly termed stall 
flutter. Since the hybrid simulation fails to predict the 
large pitching moments associated with stall, the resulting 
blade pitching of Figure 9 does not provide a representative 
picture of the stall flutter motion. 

The curves of Figure 11 illustrate the growth of blade root 
torsion with rotor lift coefficient. The hybrid results, 
compared to the measured test torsion results, were obtained 
from the theoretical elastic twist values assuming that the 
majority of the elastic twist comes from the first torsional 
mode. The trend of the hybrid results, based on full scale 
Reynolds number two-dimensional airfoil data, compares well 
with the test data except at the high thrust levels where 
large portions of the rotor are stalled. The incremental 
differences between the two curves, below stall, could be 
attributed to improper representation of the two-dimensional 
pitching moment level. A closer look at the torsional wave- 
forms for the thrust conditions above stall indicates that the 
hybrid solution does not predict the large oscillation in tor- 
sion beginning in the third quadrant. As mentioned above , 
these differences arise from an improper theoretical repre sen- 
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tat ion of the dynamic moments in the stall region. An inter- 
esting comparison is noted in Figure n between the curves of 
test data and hybrid results based on model scale Reynolds 
number airfoil data. The growth in blade root torsion ex- 
hibits a character altogether different from the test re- 
sults. This comparison indicates, as was pointed out in 
Reference 1, that even though the Reynolds number on the model 
blades is an order -of -magnitude lower than the full-scale rotor 
blades, the resulting blade loads exhibit characteristics 
similar to full scale high Reynolds number data. 

Due to the complexity of blade pressure instrumentation and 
limitations of data acquisition systems, rotor test programs 
have not recorded blade element airloads over a complete range 
of spanwise stations. Data from the BVWT 054 test is avails 
able only at the three-quarter radius station. Using the 
varied output capability of the analog system, blade element 
airloads were recorded as a function of blade azimuth and span ■ 
for each of the three advance ratios of .15, .35 and .6. 

Curves of normal force coefficient and moment coefficient span- 
wise distributions are given in Figures 12 through 19. Data 
within the reverse flow region has not been included in these 
figures. Using this data and the results from Varion traces, 
which provide similar data at each azimuth location, a summary 
of the lift stall and moment stall distributions around the 
rotor disc was made. Examples of these stall areas are illus- 
trated in Figures 19, 20 and 21. The azimuth positions for 
blade element lift and moment stall, as obtained from test 
data at the three-quarter radius, are indicated in each figure. 
The hybrid representation of the rotor tends to stall later and 
unstall earlier than the test results. 

In the unsteady aerodynamic theory, the experimentally derived 
equation for dynamic stall delay is given as: 


Dynamic stall delay = 

This equation is used in such a manner that as long as the 
time rate of change of angle of attack, oc , remains positive 
and of sufficient magnitude, stall will not occur. Thus, an 
airfoil experiencing stall delay must stall on or before the 
condition of cS: = o. This condition of ci = o has been noted 

at the three-quarter radius location on the stall regions of 
Figures 19, 20 and 21. in each case, the lift stall and moment 
stall occur prior to the condition of <=k =0. The magnitude 
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of the time rate of change of ex , therefore, is the priite 
consideration determining the azimuth location of the dynamic 
stall. As moment stall occurs on a portion of the blade, it 
develops an elastic pitching moment acting to oppose the in- 
creasing angle of attack. The rate of change of oc is, there- 
fore, reduced up to the point of lift stall and then to <k = 0. 
As cfc becomes negative, the airfoil experiences loads below 
statically predicted values. The manner in which the airfoil 
recovers from dynamic stall determines the damping character- 
istics of the airfoil. 

A measure of the damping characteristics of an oscillating 
airfoil is obtained by investigation of curves of moment co- 
efficient versus angle of attack. Reference 10 states that 
the net work done by an airfoil on the surrounding air is pro- 
portional to the area under such a curve, and is positive for 
a counterclockwise circuit and negative for a. clockwise cir- 
cuit. Thus, if the circuit encloses an area in a clockwise 
sense, it represents energy extracted from the air stream by 
the airfoil. Net energy extraction in a cycle (negative 
damping) implies that the rotor blade oscillation in which it 
occurred would tend to increase in amplitude. 

In order to visualize the effects of damping on rotor blade 
element airloads, a brief study was made using the hybrid 
simulation. In this study, the dynamic stall delay after 
ot =0, or the stall recovery, was varied by changing the V 
function in the dynamic stall delay expression. Such changes 
were made only for the portion of the cycle experiencing 
negative ot . 

The pitching moment curves of Figures 22 and 23 represent 
rotor cycles with 25%, 50%, 75% and full dynamic stall hys- 
teresis. The reduction in negative damping with decreasing 
stall hysteresis is evidenced by the reduction in area under 
the clockwise segment of the curve. The effects of improved 
damping characteristics are obvious in the curves of Figures 
24 and 25. As the damping characteristics improve, negative 
pitching moment levels after stall decrease in magnitude, re- 
sulting in reduced elastic twist and, thereby, reduced blade 
root torsion. The blade element normal force level after 
stall increases and the overall flapping excursions are re- 
duced. These effects seem to produce overall improvements in 
rotor performance and reductions in structural loads. Since 
the damping, characteristics of a rotor operating in an un- 
steady environment are a function of the choice of airfoil. 
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the question may be asked whether the characteristics of lift 
and moment stall recovery can be controlled by proper select 
tion of airfoil without seriously affecting the favorable 
dynamic lift characteristics- Experimental evidence presented 
in Reference 12 indicates that favorable dynamic lift and pitch- 
ing moment and recovery behavior can be obtained by choosing 
airfoils which do not display leading edge stall characteristics 
Since trailing edge stall airfoils do not have the high Ci max 
characteristics typical of leading edge stall airfoils, a trade- 
off must exist between improved dynamic stall recovery charac- 
teristics and reduced two-dimensional lift capability. 

To evaluate the contributions of the unsteady aerodynamics rep- 
resentation and the three-dimensional sweep effects within the 
present hybrid formulation, modifications to the present simu- 
lation were made to isolate each of these factors. Figures 26 
through 29 present a summary of the results of these modifica- 
tions for the advance ratios of .35 and .6. The quasi-static 
results were obtained by removing the sweep effects and unsteady 
aerodynamic effects from the formulation. Sweep effects were 
eliminated by removing the cosine of the sweep angle. Unsteady 
aerodynamic effects were eliminated by removing the V function 
representation of the lift and moment overshoot, and setting the 
Theororsen terms F =1.0, G * 0, eliminating the shed wake, A 
description of these terms is given in Reference 11. 
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CONCLUSIONS 


Because of its inherent computational speed and "engineer 
in the loop" on line versatility, the hybrid simulation is 
a valuable tool for studying the complex rotor environment. 
This tool can lead to a more improved representation of 
rotor unsteady aerodynamics and to the develppment of dy- 
namic design criteria for airfoils and blade elastic char- 
acteristics. 

The improvement in prediction of blade flapping in the 
stall regime over typical digital formulations may be 
attributed to the increased number of radial computations 
inherent in the hybrid simulation. A study is suggested to 
investigate the influence of radial update capability on 
prediction capability. 

The existing formulations of unsteady aerodynamics effects 
are based on oscillating airfoil data not representative of 
angle of attack oscillations experienced by the rotor. 
Studies into the effects of oscillating airfoil pitch am- 
plitude on the theoretical unsteady aerodynamic V function 
representation and the resulting effects on dynamic lift 
and moment should be made to improve existing theoretical 
methods. 

Investigation of model rotor blade root torsion results 
indicates that due to the unsteady aerodynamic environment 
of the blade, the rotor produces airloads representative of 
a full scale Reynolds number environment. 

Additional hybrid capability to represent nonuniform down- 
wash is required to predict rotor performance and blade 
airloads at the lower advance ratios ( M = .15). 
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